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摘 要 : 针对 柔性 车 间 内 AGV 最 优 蔡 换 比 问题 ， 建 立 了 员工 、AGYV 共同 搬运 的 柔性 车 间 调 度 模 型 。 该 模型 以 完工 时 
间 最 小 、 成 本 最 少 为 目标 ， 从 静态 分 析 与 动态 分 析 两 方面 求 得 最 优 替 换 比 。 静 态 分 析 部 分 主要 采用 线性 规划 完成 ， 动 
态 分 析 部 分 采用 粒子 群 算法 进行 求解 。 针 对 动态 分 析 部 分 ， 提 出 启发 式 规 则 分 配 AGV、 员 工 搬运 操作 。 通 过 算 例 得 出 
AGYV 最 优 替 换 比 的 帕 累 托 最 优 解 集 。 并 发 现 AGV 最 优 替换 比 和 AGYV 价格 有 关 。 
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Research on optimal replacement ratio of AGV in flexible workshops 


Xu Yunqin, Ye Chunming, Cao Lei 
(School of Business, University of Shanghai for Science & Technology, Shanghai 200093, China) 


Abstract: In order to solve the problem of optimal replacement ratio of AGVs in flexible workshops, this paper established a 
flexible shop scheduling model for employees and AGVs. This model took the goal of minimum completion time and cost as 
the goal and obtained the optimal replacement ratio from static analysis and dynamic analysis. The static analysis part mainly 
used the linear programming to complete, the dynamic analysis part used the particle swarm algorithm to solve. For the dynamic 
analysis part, this paper proposed heuristic rules to allocate AGVs and employee handling operations. Through experiments, this 
paper obtains the Pareto optimal solution set of optimal replacement ratio of AGV. And this paper finds that AGV optimal 
replacement ratio is related to AGYV price. 
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0 ”引言 真 来 建立 人 机 比 的 案例 研究 。 随 着 科学 技术 的 发 展 ， 越 来 越 多 
工业 机 器 人 投入 生产 使 用 ， 因 此 上 述 这 种 人 机 合作 的 模式 将 逐 

随 着 科技 发 展 ， 机 器 人 成 本 低 ， 工 作 质 量 高 ， 精 准 度 高 ， 渐 被 自动 化 蔡 代 。 
节能 减 排 等 优点 越 来 越 明 显 ， 制 造 业 或 者 服务 业 选 择 智 能 机 器 前 ， 国 内 关于 人 机 蔡 换 的 研究 主要 集中 在 从 经 济 的 角度 
人 来 代 蔡 一 部 分 甚至 全 部 人 工 ， 让 人 摆脱 简单 又 枯燥 无 味 的 操 ”分析 机 器 换 人 的 可 行 性 。 许 军 等 人 四 指出 当 机 嚣 人 成 本 逐渐 降 
作 ， 还 使 公司 获得 更 大 收益 。 最 优 人 机 替换 比 指 的 是 制造 业 投 。” 低 到 与 劳动 力 成 本 相交 的 时 候 ， 厂 商会 考虑 人 机 蔡 换 。 孙 玉 卫 


By 


资 机 器 人 代 蔡 人 工时 ， 在 可 接受 的 投资 范围 内 ， 使 得 收益 最 大 ”外 研 究 的 是 机 器 换 人 政策 补贴 以 定额 还 是 定 比 策略 补贴 给 供应 
投资 尽量 小 时 的 员工 与 智能 机 器 人 的 数量 比 。 厂商 和 应 用 企业 ， 使 得 政策 补贴 效果 最 好 。 文 献 [6,7] 研 究 的 是 
前 ,国外 关于 人 机 比 的 研究 主要 基于 机 器 不 能 自动 加 工 ， ”企业 机 器 换 人 之 后 ， 采 用 项 目 成 本 效益 去 分 析 机 器 换 人 的 经 济 


去 讨论 一 个 工人 操作 几 台 机 器 的 问题 ， 目 的 是 提高 生产 率 和 从 ”效益 。 文 献 [8] 分 析 了 汽车 厂 内 AGV 替换 人 力 前 后 的 成 本 ， 并 
动态 的 角度 ,通过 结合 调度 来 求解 最 优 人 机 比 。 文 献 [1] 通 过 工 。 分 析 影 响 成 本 的 相关 因素 。 不 难 发 现 上 述 文章 只 是 证 明了 机 器 
作 测 量 来 获取 当前 半导体 公司 人 机 比 ， 然 后 再 通过 时 间 测 量 ， 换 人 经 济 上 可 行 ， 却 忽视 研究 企业 进行 机 器 换 人 的 过 程 中 投资 
减少 不 必要 的 操作 和 取消 非 增值 活动 来 进行 工作 标准 化 从 而 优 ” 机 器 人 与 人 的 数量 比 。 
化 人 机 比 。 文 献 D2] 结 合 埃塞俄比亚 针织 服装 行业 中 不 同人 机 比 笔者 在 查阅 文献 的 过 程 发现 ， 有 学 者 研究 一 人 操作 多 机 器 
率 的 各 公司 ,去 分 析 人 机 比 的 不 同 对 裁缝 工厂 生产 效率 的 影响 。 ”时 的 人 机 比 ， 也 有 学 者 研究 人 机 蔡 换 的 经 济 性 ， 暂 时 还 未 有 学 
文献 [3] 通 过 对 约束 和 非 约束 加 工 设备 的 仿真 研究 , 提出 了 用 仿 ”者 研究 人 机 蔡 换 比 。 本 文 的 创新 点 是 将 二 者 结合 ， 使 之 更 适合 
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目前 我 国 制造 业 再 


升级 的 大 潮流 。 


在 柔性 车 间 内 ， 可 以 引进 


AGYV 替代 部 分 或 者 全 部 员工 搬运 ， 因 此 讨论 引进 多 少 AGV 与 


剩 下 多 少 员 


员工 与 AGV 的 费 | 
考虑 工件 调度 的 动态 角度 讨 
角度 进行 研究 ， 更 加 


1 


术 


员工 搬运 。 
柔性 
工 年 的 工资 收入 为 8 元 , 辞退 
搬运 速度 v, 件 /min。 现 有 平均 月 订单 数 m 每 天 工作 8 小时, 二 


单 的 准 
使 人 、 机 利 


工 有 一 定 的 现实 意义 。 


本 文 研究 柔性 车 间 AGYV 最 优 蔡 换 比 问题 时 , 以 完工 时 间 、 


确 预 测 。 如 果 订 单 预 测 过 
j 率 变 低 。 


1 柔性 车 间 内 AGYV 最 优 蔡 换 比 建 模 


为 目标 ， 从 不 考虑 工件 调度 的 静态 角度 和 
论 最 优 蔡 换 比 。 
符合 现实 情况 ， 可 以 作为 企业 进行 车 
几 蔡 换 的 参考 依据 。 在 企业 做 人 机 蔡 换 决策 
寸 少 , 造成 违约 成 本 , 预测 过 多 ， 


从 微观 和 宏观 两 个 
间 人 
时 ， 关 键 在 于 对 订 


某 公 司 想 在 柔性 生产 车 间 内 引进 AGV 代替 部 分 或 者 全 部 


由 于 流动 资金 有 限 ， 和 希望 投资 后 成 本 低 于 投资 前 。 
车 间 原 有 X 人 负责 搬运 ， 引 进 AGV 后 剩 下 x 个 工人 。 


Pal 


inaXiv 合 作 期 十 


| 
等 : 柔性 车 间 内 AGV 最 优 替换 比 研究 


前 提 下 , 使 成 本 最 小 的 人 机 组 合 。 那 么 , 最 优 蔡 代 比 = 员工 数量 : 
AGYV 数量 =0: 3， 且 较 投资 前 ， 工 作 时 间 缩 短 3.61 天 ， 成 本 减 
少 149.8 万 。 
[3 0 18.06 120.2 | 
4 0 13.54 151.1 
5 0 10.83 182 
6 0 9.03 212.9 
[AGV 人 天 数 成 本 ]=|7 0 7.74 243.8 
3 1 17.11 137.7 
4 113 168.6 
5 1 10.48 199.5 
[6 1 8.78 230.4 | 
1.3 单价 和 工资 敏感 性 分 析 
由 AGYV 价格 分 析 可 知 ，AGYV 价格 跨度 很 大 ， 有 些 低 至 5 
万 ， 有 些 高 达 100 万 ， 因 此 讨论 不 同 价格 的 AGV 对 最 佳 蔡 换 
比 的 影响 是 有 现实 意义 。 由 下 图 可 知 ， 当 AGYV 单价 在 10 到 50 
万 时 ， 最 优 人 机 蔡 换 比 为 0:3; 当 AGYV 单价 在 60 到 80 万 时 ， 


一 个 员工 需要 ci 元 。 工人 的 平均 


最 优 人 机 蔡 换 比 为 2:3; 当 


an 元 。 


(D) 


班 制 , 每 人 每 月 排 班 工作 22 天 。 购买 AGV 需要 Psov 元 / 台 , 使 
用 寿命 n 年 ,， 共 y 台 。AGYV 的 平均 搬运 速度 vaov 件 /min， 维 修 
费用 wei 元 /年 / 台 。AGYV 投入 使 用 前 的 安置 费用 
一 | m 
目标 函数 : hi =min( 
2。yeyAGv 十 XeV 人 
户 =min(CAcv 二 CN) 
约束 条 件 : 
CovtC, <Xxg, Xn 
m m 
2XxyxyaGcv 十 XXxVyN XXvVN 
Cov = yx paov +ant+wei 
C=xxg, Xn+(X—x)xci 
xeN',yeN 
x<X 
1.2 问题 求解 


员工 需要 付 一 次 性 辞退 费 0.5 


投资 前 ,15 个 员工 负责 搬运 。 员 
万 。 


即 160 件 / 天 。 
投资 后 ，AGYV 单价 30 万 每 台 ，AGYV 搬运 速度 1 件 /min， 


每 天 不 轮班 ， 即 960 件 /天 。 


资 每 年 6 万 , 辞退 一 个 
员工 搬运 速度 为 1 件 /3min， 


一 次 性 总 安装 费用 


20 万 ，AGYV 维 


修 费 每 台 每 年 0.3 万 。 且 投资 前 后 ， 和 柔和 


§ 


为 该 柔性 车 间 引 进 AGV 前 ， 
工 速度 是 处 于 一 个 平衡 状态 。 当 车 间 只 3 
产品 速度 会 受到 加 工 设备 的 约束 ， 因 此 选择 能 完成 该 订单 


2000 件 。 

通过 MATLAB ji 
* AGYV 的 数量 、 员 
资 前 


运算 得 到 以 下 非 劣 解 。 
的 数量 、 需 要 


， 并 提高 整体 搬运 速度 ， 


车 间 平均 


月 订单 件 


其 中 每 一 列 依次 表 


从 而 帮助 企业 者 


作 的 天 数 、 成 本 。 那 么 投 
和 [0 15 21.67 270] 由 结果 可 以 看 出 ， 引 进 AGV 可 以 减少 成 
得 更 大 的 利润 。 
各 加 工 设 备 加 工 速度 与 搬运 的 员 
进 AGV 时 ， 整 体 的 


性 车 间 不 应 引进 AGV 。 


AGYV 单价 大 于 等 于 90 万 时 ， 该 柔 


所 不 同 ， 


由 员工 工资 具体 行业 工资 不 同 , 且 随 着 时 间 变 化 
因此 有 必要 研究 员工 工资 对 最 优 人 机 蔡 换 比 的 影响 。 


可 以 看 出 工人 年 工资 在 6 万 到 10 万 之 间 ， 该 柔性 车 
最 优 蔡 换 比 不 变 ， 均 为 0: 3。 


10 20 30 40 50 60 70 80 90 100 


由 图 1、2 


数量 
和 上 人 -AGV 
12 
10 
8 
6 
4 一 一 一 一 二 
-4 4 刁 之 22 
, 了 一 陡 一 隆 -了 dn 
“时 中 嘻 叶 ET 
学 
(从 


间 的 人 机 


图 1 AGYV 单价 /万 与 最 优 人 机 蔡 代 比 (工资 -6 万 /年 ) 


3 

4 

3 人 机 比 =0: 3 

2 

1 

0 

6 6.5 7 7.5 8 853 9 9.5 10 

图 2 工人 年 工资 /万 与 最 优 替 代 比 (AGV=30 万 / 台 ) 


2 ”动态 分 析 一 一 AGV 最 优 替 换 比 


2.1 柔性 车 间 内 AGV 最 优 替 换 比 问题 描 
2.1.1 问题 描述 

柔性 车 间 AGV 最 优 蔡 换 比 问题 考虑 原 有 X 个 员 
运作 业 ， 现 因 人 工 成 本 过 高 的 问题 准备 投资 AGV 来 负 


述 与 建 模 
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搬运 作业 并 解雇 部 分 员 


。 投资 AGV 后 员工 与 AGYV 的 总 成 本 
指 的 是 从 AGV 投入 使 用 到 报废 整个 周期 内 , 员工 相关 费用 ( 含 
工资 、 辞 退 费用 ) 及 AGV 相关 费用 〈 含 AGV 产品 成 本 、 相 关 
安装 费 、 维 修 费 )。 投 资 AGV 前 总 成 本 指 的 在 相同 时 间 内 ， 员 
工 相 关 费 用 ( 含 工 资 、 辞 退 费 用 )。7 工件 在 mm 台 设 备 上 加 工 。 
x 个 员工 与 y 台 AGV 共同 负责 工件 在 设备 间 的 运输 。 且 
每 个 员工 运输 和 AGYV 运输 是 相互 独立 的 。 每 个 零件 都 及 道 
不 同 的 工序 %<m , 每 个 工序 可 在 多 台 设 备 上 加 工 。 有 如 下 假设 : 
a)AGYV 的 搬运 路 线 是 固定 的 ,AGYV 搬运 不 延 时 ,不 同 AGV 
搬运 互 不 干扰 。 
b) AGYV 工作 内 容 是 : 从 上 一 


道 工序 加 工 设备 的 缓存 区 立 


即 搬运 至 下 一 道 工序 的 缓存 区 

c) 员工 的 搬运 路 线 是 固定 的 ， 员 工 搬运 不 延 时 ， 不 同 员工 
搬运 互 不 干扰 。 

d) 员工 搬运 工作 内 容 是 : 从 上 一 道 工序 加 工 设 备 的 缓存 区 
立即 搬运 至 下 一 道 工序 的 缓存 区 


e) 零 时 刻 ， 所 有 设备 和 机 器 都 可 用 

f) 工件 从 加 工 设备 的 缓存 区 到 开始 加 工 的 时 间 是 瞬时 完成 ， 
消耗 时 间 忽略 不 计 

g) 每 台 设 备 一 次 只 和 
工时 不 中 断 。 

h) 设备 加 工 顺序 约束 ， 加工 设 备 按 先 到 先 加 工 (FCFS) 规 


能 加 工 一 个 零件 , 每 个 零件 在 设备 上 加 


则 。 


i) 不 同 工 件 先后 无 约束 , 同 种 工件 必须 先 完成 上 一 道 工序 ， 
然后 再 进行 下 一 道 工序 。 

j) 所 有 待 加 工 的 工件 放 在 装载 台 M0， 加 工 好 的 工件 放 在 
纯 载 台 Mm+l 处 。 


2.1.2 模型 建立 
Cov 表示 AGV 相关 费用 ( 含 AGV 产品 成 本 、 相 关 安 装 
费 、 维 修 费 ) 


CAX 表示 员工 相关 费用 〈 含 工资 、 

n 种 工件 ， 件 i{1,2,3,.…n} 

m 台 设 备 ， 装载 台 Mo ， 
j:{0,1,2,3,...m,m+1} 

X 个 员工 ， 投 资 前 员工 数量 

x 个 员工 ， 投 资 后 员工 Human:{1,2,3,...x} 

y 台 AGV， 投 资 后 搬运 机 器 人 R:{1,2,3,.…y} 

每 个 工件 i1 有 nn 道 工序 , 工序 k:{1,2,3,…n,} ， 即 工件 i 的 第 
k 道 加 序 : 0O. 

各 工序 0 开始 加 工时 刻 ，# 工序 Oi 完成 加 工时 刻 。 

Pix 工序 Oi 的 加 工时 间 。 

ht 工序 0,4_ 完成 后 ,员工 开始 搬运 工件 至 工序 
工 设备 缓存 区 的 时 刻 。 

Mix 员工 完成 搬运 时 刻 。 

rt 工序 Oi 完成 后 ，AGV 开始 搬运 工件 至 工序 0 所 在 
加 工 设 备 缓存 区 的 时 刻 。 


辞退 费用 ) 


色 载 台 M,， 设 备 


Wy 


Ox 所 在 加 


rt 表示 AGV 完成 搬运 时 刻 。 
六 表 示 员 工 从 设备 方 搬运 工件 至 设备 户 所 需要 的 时 间 


j7,。 表 示 员 工 从 设备 广 搬 运 空 载 至 设备 户 所 需要 的 时 间 


rTi.; 表示 AGV 从 设备 ji 搬运 工 


件 至 设备 户 所 需要 的 时 间 


rT ,表示 AGV 从 设备 方 搬运 空 载 至 设备 户 所 需要 的 时 间 


标 函 数 =min(maxti,) 
万 =min(CAcv 十 CN) 
约束 条 件 : 
a) 工件 从 装载 台 搬 运 出 来 进行 第 一 道 工 序 的 加 工 , 且 工序 
Oi 1 在 加 工 设备 j 上 加 工 。 
Sp, be 2) 
yp, #1 十] Sh (3) 
b) 工序 O; 在 加 工 设备 j 上 加 工 。 
和 Pe, ) Si (4) 
c) 搬运 完成 之 后 才能 进行 加 工 。 
Ogre < (5) 
0ghr. St (6) 
d) 工序 Op 加 工 完成 ， 立 即 搬运 。 
fi < ri (7) 
ti < ht (8) 
e) 工 件 前 后 工序 搬运 时 间 约 束 : AGV 或 者 员工 从 工序 O; 4 
所 在 加 工 设备 搬运 至 工序 O; ,1 所 在 加 工 设备 。 
fin = -> pj Np *77; j 十 7 HL (9) 
Ha 二 > San, * nj #7 j ,+ ht KH (10) 


f) 相 邻 两 道 搬 运 任务 时 间 约 束 : 一 个 AGV (员工) 完成 上 
道 搬运 任务 ， 行 驶 到 下 道 搬 运 任务 的 初始 加 工 设 备 ， 再 将 该 工 
件 搬运 至 加 工 设备 上 。 


汪 T 


DD ED 六 僻 帮 1 
所 名 名 名 (GD 
*, koa rT Sn w+MI( 1 pp 
需 需 部 总 
2 jl. 
hi + 2 Dn 于 区- 大 AhT) ) + 
所 所 名 所 (12) 


TEN 二 MG 一 Ne) 
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g) 同一 加 工 设 备 加 工 的 工序 约束 : Onn 、 Op 先后 在 同一 


加 工 设备 加 工 。 
而 让 + Pip Sito +M(l— Zi nbih) (13) 
h) 一 道 工序 只 在 加 工 设备 加 工 一 次 。 
2 xuj=1 (14) 
六 
i) 工序 Oo 在 加 工 设备 j 上 加 工 ， 搬 运 至 卸载 台 。 
帮 ， + > Pin we +77) fin (15) 
六 
i Sp, sx j+hT,,, <t (16) 
六 
j) 其 他 。 
M=+o0 (17) 
1 若 工 序 O,, 在 设备 j 上 加 工 
Nip,j = 0 此 余 (18) 
1 若 工 序 O,; 和 0O, ,由 同一 AGV 搬 运 
WYip bh = 0 其 余 (19) 
jw 人 考 工 序 0 和 0 外 同一 员 T 扳 运 。 20 
交办 0 其 余 ) 
1 若 工序 0 和 O, 在 同一 设备 加 工 
Zi bb = 0 3 Se (21) 
余 


2.2 构建 帕 累 托 最 优 解 集 
如 果 调 度 方案 x 为 帕 累 托 最 优 解 ， 表 示 方 案 x 不 被 可 行 调 


度 集 合 中 任何 一 个 调度 方案 所 支配 ， 
函数 值 必 须要 削弱 至 少 一 个 其 他 目标 函数 值 踢 。 方 


即 如 果 要 改进 x 对 应 的 任 


目标 函数 值 不 比方 案 y 


何 一 个 目标 
案 x 支配 方案 y 表示 方案 x 对 应 的 所 有 
对 应 的 所 有 目标 函数 值 差 ， 日 


至 少 存在 一 个 目标 函数 值 优 于 方 


案 y 所 对 应 的 目标 函数 值 。 多 目标 调度 问题 的 最 优 解 不 是 唯一 


的 ， 而 是 
本 文 结合 双 元 和 


90。 首先 取出 方案 集合 A 中 第 一 个 方案 x， 剩 下 的 方案 组 成 


多 个 帕 累 托 最 优 解 组 成 的 帕 累 托 最 优 解 集合 。 


负 标 赛 与 动态 淘汰 制 来 构建 帕 累 托 最 优 解 生 


mr lt 


浪 


合 B。 将 方案 x 与 外 
合 B 中 不 存在 支配 
| 除 ， 


蚌 合 B 中 的 方案 一 一 比较 。 如 果 剩 余 方案 集 
x 的 方案 , 那么 把 B 中 被 x 方案 支配 的 方案 


且 将 方案 x 放 入 帕 累 托 最 优 解 集合 C 中 。 如 果 剩 余 方 案 


Rin 


案 B 中 任何 一 个 方 


中 存在 方案 支配 x 方 案 , 那么 将 x 方案 从 方案 集合 A 中 
且 删 除 剩余 方案 B 中 被 x 支配 的 方案 。 


如 果 x 与 剩余 方 
案 都 互 不 支配 ， 那 么 从 方案 集合 A 中 删除 


将 x 放 入 集合 C 中 。 然 后 取出 方案 集合 A 中 第 二 个 方案 ， 


以 出 


上 类推， 直到 方案 集合 A 为 空 。 


采用 该 方法 构造 帕 累 托 最 优 


解 集 可 以 在 每 个 循环 中 可 以 实时 淘汰 被 支配 解 ， 从 而 达到 降低 


复杂 度 的 
2.3 


目的 。 


粒子 群 算法 设计 
柔性 车 间 调 度 问题 往往 是 NP 难 问题 , 而 含有 AGYV 的 柔性 
车 间 调 度 问 题 在 柔性 车 间 调 度 问 题 基础 上 考虑 了 AGYV 分 配 ， 


使 得 问题 更 难 求解 。 


含有 AGYV 的 柔性 车 间 调 度 问 题 往往 很 难 


通过 精确 的 数学 方法 求 得 最 优 解 ， 而 是 一 般 使 用 启发 式 智能 算 


法 得 到 相对 最 优 
术 , 具有 较 高 的 


行 性 ， 


。 粕 子 群 算法 是 一 种 基于 种 群 的 智能 优化 技 


次 运行 可 以 搜索 解 空间 中 的 多 个 解 ， 


因而 非常 适用 于 求解 柔性 车 间 调 度 问题 。 


2.3.1 算法 编码 
本 文 用 
3L, 1 


的 。 与 让 


个 粒子 位 置 代 表 
Tr (transport sequence) 向 
(Crobot-human ) 向 量 构 成 ， 日 


这 


个 候选 


三 
EE 


EM 


周 度 方 案 ， 其 中 维 
(machine) 向 量 、RH 


个 向 


量 每 


维 上 是 相互 对 应 


1 面 两 章 不 同 的 RH 向 量 表示 含 AGV、 员 工 的 搬运 向 量 。 


假设 员工 x=1，AGYV 数量 y=1 时 ， 下 表 为 一 个 编码 示例 。 那 么 


RH 向 


三 | 
星 T， 


1 表示 员工 搬运 ，2 表示 月 


表 1 


粒子 编码 向 量 和 矩阵 


有 AGYV 搬运 。 


维 数 


10 


1. 12 


13 14 15 16 


1 
Tr 向 量 1 
M 向 量 1 
RH 向 量 1 


PD bb bb 
一 mw w|w 
ID 上 上 三 
天 一 iu 
ND w w|c 


王 DI 


9 
4 
1 
1 


DD mm 一 |o 


1 
5 
2 


3 
4 


[> 


1 


» 
a 
a 

iD 只 


2.3.2 种 群 初 始 化 


根据 上 文 编码 方式 ， 
献 [11] 初 始 种 群 定位 法 思想 ， 


站 乡 


E 随 机 生成 搬运 序列 。 然 后 借鉴 文 


选择 操作 工序 对 应 加 工时 间 最 短 


加 工 设备 ， 得 到 M 向 量 。 


的 加 工 设 备 。 最 后 利用 启发 式 规 则 生成 RH 向 量 。 


a) 算法 实现 随机 生成 多 个 搬运 序列 即 Tr 向 量 。 


具体 步骤 如 


b ) 为 每 个 搬运 序列 所 对 应 的 加 


c) 为 每 个 搬运 序列 的 搬运 了 


[ 序 


序 分 配 加 工时 间 最 短 的 


启发 式 分 配 AGV、 员 工 搬 


运 ， 得 到 RH 向 量 。 启 发 式 规则 :选择 最 早 能 够 完成 该 搬运 工序 


的 人 或 AGV 进行 搬运 即 min{hti ,rt 1} 。 
备 加工 完 成 后 ， 需 要 搬运 至 加 工 设备 j 加 工 0;; ， 分 配 最 早 


能 够 到 达 的 员工 或 者 AGV 进行 搬运 。 
2.3.3 粒子 位 置 更 新 


本 文通 过 同时 结合 文 圳 
出 的 RPX 交叉 ， 提 出 了 粒 了 
想 可 知 ， 粒 子 的 位 


大 [12] 提 


当 工序 0 在 加 工 设 


上 的 更 新 方法 与 文献 [13] 提 


群 


位 


Sp 


其 中 :%@,， 06, Cc， ,号 
k 次 迭代 时 单个 粒子 的 历史 最 优 位 置 ， g3; 表示 第 k 次 迭代 时 


全 局 最 优 位 置 


与 基 


更 新 方法 。 
是 “自身 “ 认 知 “社会 "三 部 分 互相 影响 作 
的 结果 。 因 此 定义 粒子 位 置 更 新 公式 如 下 : 

Xi! = C1, ® hor ® al ® fxy, pB!), eB!) 


车 本 粒 子 群 思 


(22) 


本 粒子 群 算法 一 致 ，pB; 表示 第 


身 部 分 


FF=w® | 


f(x) 


Kk 
在 
Es 


rand<w 


身 部 分 表示 粒子 对 上 


次 


a 


，R 向量 跟 
a) 交 换 。 随 机 选择 两 个 不 同 


随 Tr 


otherwise 


向 量变 换 。 


(23) 


身 的 思考 ， f(x’) 表示 粒子 
当 rand 产生 的 随机 数 小 于 。 则 执行 以 下 操作 ， 否 则 保持 

第 一 步 变换 Tr 向 量 : 从 以 下 三 个 操作 随机 选择 一 个 执行 。 
且 在 执行 操作 时 ，M 


件 的 


序 ， 


交换 两 个 基因 人 


小 
3 
局 
并 
稿 
人 
ylll 
Da 
i 
Es 
现 
(ULD 
到 


图 3 Tr 向量 交换 
b) 插 入 。 随 机 选择 两 个 位 置 a 和 b, 将 b 基因 插入 都 的 前 
面 ，a 到 b 之 间 的 基因 后 移 一 位 。 操 作 示 意图 如 图 4 所 示 。 


pa 
ys 
二 
只 


图 4 Tr 向 量 插入 
0O) 反 序 。 随 机 选择 两 个 位 置 a 和 b, 将 a, b 之 间 (包括 a， 
b) 之 间 的 基因 序列 反 转 ,其 余 片 段 不 变 。 操作 示意 图 如 下 图 所 


tI 


图 5 Tr 向量 反 序 

第 二 步 变换 M 向 量 。 随 机 选取 一 个 位 置 a， 随 机 从 该 工 
可 行 加 工 加 工 设备 集 选 取 另 一 个 加 工 设备 蔡 换 a。 如 果 第 一 步 
随机 选择 的 是 交换 操作 ， 则 M 向 量 及 M 向 量变 换 操作 示意 
如 图 6 所 示 。 


I 


= 


图 6 M 向 量变 换 
第 三 步 变 换 R 向 量 .随机 选取 一 个 位 置 a, 随 机 从 搬运 AGV 
兰 合 中 选取 另 一 个 AGV 替换 a。 操 作 示 意图 如 图 7 所 示 。 


ye 


2) 认 知 部 分 


SC PDB) rand<or 


0 = 1@S0 pB;)= | (24) 
i otherwise 


认 知 部 分 表示 粒子 对 历史 位 置 的 思考 。 当 rand 产生 的 随机 
数 小 于 cz 则 执行 POX 交叉 操作 ， 和 否则 保持 上 一 步 自身 部 分 操 
作 结果 不 变 。 随 机 选择 工件 a， 经 POX 交叉 所 得 两 个 子 代 ， 随 
机 选择 一 个 作为 ， 操 作 示意 图 如 图 8 所 示 。 

3) 社会 部 分 


NG 2B1) rand< c,, 
CG 


Wan @ ho a0)-| (25) 


otherwise 


社会 部 分 表示 粒子 对 全 局 位 置 的 思考 ， 同 时 展现 了 粒子 之 
闻 的 协作 。 当 rand 产生 的 随机 数 小 于 6 则 执行 RPX 交叉 操 
作 ， 和 否则 保持 上 一 步 自身 部 分 操作 结果 不 变 。 随 机 产生 长 度 与 
搬运 序列 相同 的 Rand 序列 , 计算 自 适应 调整 概率 (其 中 t 为 当 


徐 云 至 ， 等 : 末 仁 车 为 AGV 人 尝 机 比 研 究 


前 迭代 次 数 ，T 为 总 次 数 ，, 为 Rand 序列 最 大 值 ，,, 为 最 


小 值 ), 经 RPX 交叉 得 子 代 1; ， 操 作 示 意图 如 图 9 所 示 《〈 假 设 


伍 0.7)。 自 适应 调整 概率 表达 式 如 下 : 


9 ”RPX 交叉 


2.3.4 算法 求解 步骤 


算法 步骤 如 图 10 算法 步 又 图 所 示 。 
的 循环 伪 代 码 如 下 : 
初始 化 员工 数量 为 -1，AGV 数量 为 9 
while( 投 资 后 AGV 总 成 本 与 人 工 总 成 本 <= 


员工 数量 = 员工 数量 +1; 


while( 投 资 后 AGV 总 成 本 与 人 工 总 成 本 <= 投 资 前 人 工 总 成 本 ) 


AGV 数量 =AGV 数量 +1; 


投资 后 总 成 本 =AGV 总 成 本 安装 费 
人 工 总 成 (辞退 费 、 员 工 工资 ); 
求 出 此 时 最 短 完工 时 间 ; 
end 
初始 化 AGV 数量 为 6; 
end 
2.4 问题 求解 
2.4.1 算法 寻 优 能 力 验 证 


实验 仿真 环境 为 Windows7 64 操作 系统 ， 处 理 器 了 


(26) 


其 中 求解 AGV 替换 比 
投资 前 人 工 总 成 本 ) 


、 维 修 费 、AGV 总 价 ) + 


E 频 2.5 


GHz， 内 存 4 GB， 采 用 MATLAB R2014 编译 软件 ， 每 个 算 侦 
运行 20 次 取 最 优 值 。 在 粒子 群 算法 中 有 五 个 超 参 数 需要 选择 ， 


分 别 是 粒子 群 种 群 规模 ， 和 迭代 次 数 ， 惯 怕 


常量 以 及 两 个 学 习 


大 
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子 。 本 文中 的 学 习 因 子 和 惯性 常量 的 选择 主 


[9,12,14~18]， 选 取 学 习 因 子 c、6 为 0.9， 惯 性 常量 中 取 0.9。 
粒子 群 种群 规 模 和 从 代 次 数 


合理 值 ， 最 后 选取 粒子 群 种 群 规 模 为 500， 迭 代 次 数 为 100。 车 间 刘 
FE 粒子 群 算法 的 寻 优 能 力 ， 从 Bilge 等 人 09 提 出 的 
选取 部 分 进行 计算 : 4 台 加 工 设 备 、 装 载 台 、 务 载 
局 ， 并 且 考 虑 AGV 的 搬运 时 间 和 空 
。 算 例 “Ex31” 表 示 ,3 表示 Jobset3， 


为 了 验 详 
40 个 算 例 
台 , 与 AGV 构成 了 四 
载 时 间 , 求 出 最 短 加 工时 间 


种 布 


要 是 根据 计算 机 的 处 


要 参照 参考 文献 完 
所 得 到 
速度 选取 。 结果 。 


| 的 结果 。FDE、FMAS 分 别 表示 文 击 


成 搬运 至 站 载 台 的 情况 下 ， 采 用 遗传 算法 与 启发 式 混 合算 法 
[17, 18] 提 出 算法 的 


PSO 和 FPSO 表示 使 用 本 文 改 进 粒 子 群 算法 ， 
蛋 度 和 和 柔性 车 间 调 度 的 最 短 加 


本 文 提出 的 改进 粒子 群 算法 有 较 好 的 寻 优 能 力 : 运 | 


1 表示 layoutl。MAS 表示 文献 [15] 提 出 算法 的 


GAHA-3 表示 文献 [16] 


算 例 


EX11 
EX12 
EX13 
EX14 
EX21 
EX22 
EX23 
EX24 
EX31 
EX32 
EX33 
EX34 


MAS-2 
[106] 
R=2 
130 


(加 工 
短 


(最 早 


r: 随机 生成 掀 
运 序列 
MA: 选 择机 床 


RH 分 配 AGV、 人 


疆 电 


结果 。 GAHA-2、 
是 出 的 2 个 和 3 个 AGYV 在 考虑 工 从 


时 间 最 
) 


完成 ) 


初始 化 


加 工 


TT 


工时 间 。 通 过 实验 结果 可 知 


分 别 计算 


本 文 算法 


计算 车 间 调 度 问 题 时 ，66.67% 的 解 优 于 MAS (R=2) 与 GAHA 


(R=2) 的 最 优 解 ，100% 的 解 优 于 GAHA (R=3) 外 
J 


算 和 柔性 车 


9 解 。 在 计 


间 调 度 问 题 时 ，50% 的 解 优 于 FDE (R=2) 与 FMAS 


(R=2) 的 最 优 解 。 
工时 间 ， 完 工时 间 越 长 ， 增 加 AGYV 更 


1234232141323413 
1234152315463666 
1222121112121212 


SB = min(x) 


1214232341323413 
1231152315463666 
1222121212121212 


种 群 


1234231421313423 
4131522431464666 
i212122112121212 


1234232141321343 
4231122241466466 
T212121212122122 


认 知 


解 的 评价 


自身 


1214232341323413 
1231152315463666 
1222121212121212 


1213423241323413 


1213244321323413 
4133114331464666 
i221121212121212 


1214233421334123 
1131553135436666 


1222122112112112 


RPX 


产生 新 的 种 群 


迭代 次 数 =t 


POX 


最 大 迭代 次 数 


1231152315463666 
1222121212121212 
1212324341323413 
1231152315463666 
1222121212121212 


x 个 员工 


y 台 AGV 
最 优 解 


同时 不 难 发 现 ， 增 力 


r= 
XX o 


1 AGYV 能 够 缩短 整个 完 


其 中 ， 如 下 图 所 示 为 R=2 时 ，FDE、FMAS、FPSO 三 种 算法 求 


X=X+l 
解 的 评价 
y=y+1 否 
双 元 锦标 赛 帕 累 托 最 优 
动态 淘汰 制 解 集 换 比 
图 10 算法 步骤 图 
表 2 算 例 试验 结果 
GAHA-2 Pso-2 FDE FMAS fpso-2 GAHA-3 Pso- 解 各 算 例 的 最 小 加 二 时 间 » 

[27] 第 三 章 [107] [108] 第 [27] 第 150 
R=2 R=2 R=2 R=2 R=2 R=3 R=3 
116 116 96 111 99 96 90 130 
91 92 82 87 80 86 79 
100 98 84 91 83 90 90 110 
140 139 103 128 108 110 110 
121 125 100 128 101 104 104 90 
88 92 76 88 74 84 81 
95 94 86 102 87 86 86 70 
151 162 106 131 114 118 116 
134 98 99 114 96 105 87 
101 79 85 99 76 89 74 
107 81 86 102 79 91 74 
166 133 110 128 i27 130 100 


图 11 


R=2 时 算法 Makespan 折线 图 
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2.4.2 示例 求解 


为 了 探究 AGV 最 优 蔡 换 比 ， 求 解 如 下 示例 。 


投资 前 ,15 个 员工 负责 


举例 第 二 行 第 三 列 数值 18= 算 
从 设备 1 出 发 到 达 设 备 2 需要 


36 
54 
39 
36 


18 
$1 


投资 后 , AGYV 搬运 时 间 久 
一 列 分 别 表示 装载 台 、 介 载 台 。 举 例 第 二 行 第 三 列 数值 6= 第 三 


发 到 达 设 备 2 需要 6 单 


求解 得 到 帕 累 托 最 优 解 入 


12 


0 24 36 
24 0 18 


0 
57 
24 
21 
27 
Sl 
72 
39 


63 3 
0 24 36 


运 。 员工 工资 每 年 6 万 。 搬运 时 间 
和 矩阵 如 下 ， 其 中 第 一 列 和 最 后 一 列 分 别 表 示 装 载 台 、 彼 载 台 。 


51 
27 
54 39 36 


三 列 行 


45 
51 


66 27 63 


42 
57 
69 


第 二 行 数 值 18， 表 示 员 工 
18 单位 时 间 。 


0 33 54 
33 0 51 
54 51 0 75 
60 75 0| 


E 阵 如 下 所 示 ， 甚 中 第 一 列 和 最 后 


车 间 加 工人 各 


列 第 二 行 数值 6, 表示 AGV 从 设备 1 
位 时 间 。AGYV 单价 30 万 每 台 , 一 次 性 总 安装 费用 
维修 费 每 台 每 年 0.3 万 , 辞退 一 个 员 了 
万 。 且 投资 前 后 ， 柔 性 


17 20 


19 0 
8 
和 
9 4 
17 
24 6 21 
13 
21 
12 


18 13 
19 8 7 


12 0 10 


9 


15 


12 22 9 21 
16 4 5 
17 24 
l]2 5 4 9 6 


12 


5 10 0 5 23 


15 


3 0 2 


9 12 23 2 0 


18 
13 4 17 
10 18 9 

18 13 


工 数量 ， 最 短 完 工时 间 ， 总 成 本 ]。 
工 规 模 上 进行 替换 ， 即 投资 后 总 成 本 小 了 
因此 AGYV 最 优 替 代 比 = 员工 数量 : AGV 数量 =1: 7， 


12 
15 
17 


E 阵 如 下 ， 


14 0 


19 
23 


12 22 9 21 


12 


20 万 , AGV 


[需要 付 一 次 性 辞退 费 0.5 
E 务 为 15X10[96]。 


20 25 0 |] 
14 19 8 
13 21 12 
13 10 18 
21 4 18 13 
18 17 9 12 
15 17 22 
19 23 9 
11 18 21 
11 0 17 20 
18 17 0 25 
20 25 0 | 


12 


14 


且 和 矩阵 


[AGYV 数量 ， 员 


天 


为 该 柔性 车 间 是 在 现 有 加 


F 等 于 投资 前 总 成 本 。 


较 投 资 


前 [0 15 380 270] ， 加 工 相 同 任务 15X10[96] 最 短 加 工时 间 缩 
短 105 单位 时 间 ， 成 本 减少 8.7 万。 


帕 累 托 最 优 解 = 


i 


OO 一 Om 人 人 Or- 


Lo 


1603 58.4 | 


816 
561 
439 
355 
307 
285 
275 


1234 75.9 


744 
539 
353 


89.3 
120.2 
151.1 
182 
212.9 
243.8 
261.3 


106.8 
137.7 
199.5 | 


换 比 , 没有 考虑 了 
有 结合 具体 工厂 实 


2.4.3 单价 和 工资 敏感 性 分 析 


AGYV 单价 跨度 较 大 ，| 


轻 负载 AGV 单价 


第 二 章 搬运 机 器 人 市 场 分 析 可 知 ， 
5-8 万 , 中 高 端 AGV 单价 20-50 万 , 甚至 100 


万 。 因 此 有 必要 研 


下 表 可 以 看 出 当 AGYV 单价 小 于 50 万 元 时 ,车 间 可 以 采用 
去 替代 人 工 ， 否 则 采用 全 是 人 


究 AGV 单价 对 最 优 人 机 蔡 换 比 的 影 1 


= 


。 | 
AGV 
司 时 可 以 看 出 AGV 单 


更 好 。 


价 对 蔡 换 比 的 影响 较 大 ， 且 AGYV 的 单价 越 高 ， 车 间 应 该 减少 


投资 AGV 数量 。 


表 3 AGYV 单价 (万 ) 与 替代 比 (工资 =6 万 /年 ) 


AGV 单 价 AGV 人 完工 时 间 成 本 
10 20 0 187 245.5 
15 15 0 201 266 
20 11 0 219 257.4 
25 9 0 248 260.6 
27 8 0 254 250.7 
30 7 275 261.3 
35 6 0 307 242.9 
40 5 j 353 249.5 
45 5 0 355 257 
50 0 15 380 270 
55 0 15 380 270 
员 资 具 体 行 业 工 资 不 同 ， 且 随 着 时 间 变 化 有 所 不 同 ， 
因此 有 必要 研究 员工 工资 对 最 优 人 机 替换 比 的 影响 。 由 下 表 可 
以 看 出 工人 工资 对 替换 比 的 影响 较 小 。 
表 4 工人 工资 (万 /年 ) 与 替代 比 (AGV=30 万 / 台 ) 
工资 AGV 人 完工 时 间 成 本 
6 了 1 275 259.8 
6.5 T 1 275 259.8 
8 了 1 275 259.8 
7.5 了 1 275 259.8 
8 7 1 275 259.8 
8.5 7 1 275 259.8 
9 7 0 285 243.8 
9.5 7 0 285 243.8 
10 了 0 285 243.8 
3 ”结束 语 
本 文 创新 性 地 提出 了 柔性 车 间 AGV 替换 比 模型 ， 并 从 静 


动态 两 个 角度 


去 分 析 最 优 人 机 比 。 更 贴 合 实际 生产 加 工 情 


具有 很 强 的 参考 及 应 用 价值 。 


本 文 在 从 静态 4 


加 工 设 备 的 约束 ， 

人 机 比 即 为 最 优 人 
工时 间 最 短 角度 出 
投资 前 的 基础 上 , 


度 分 析 时 ， 主 要 从 成 本 角度 出 发 ， 考 虑 到 
在 满足 加 工 任务 的 基础 上 ， 对 应 成 本 最 小 的 
机 蔡 换 比 。 在 从 动态 角度 分 析 时 ， 主 要 从 加 
发 ， 考 虑 加 工 设 备 的 约束 ， 在 满足 成 本 低 于 
应 加 工时 间 最 短 的 人 机 比 即 为 最 优 人 机 比 。 


从 动态 角度 分 析 时 
行 的 生产 调度 ， 其 
本 文采 用 针对 性 的 
员工 与 AGV 搬运 
锦标 赛 与 动态 淘汰 


度 。 


， 本 文 创新 性 地 采用 3xL 维 向 量 表示 一 种 可 
中 RH 向 量 表示 人 、AGYV 的 搬运 情况 。 同 时 
粒子 群 更 新 规则 ， 同 时 结合 启发 式 规则 分 配 
， 从 而 算法 得 到 更 优 的 结果 。 最 后 结合 双 元 
制 来 构建 帕 累 托 最 优 解 集 ， 减 少 了 算法 复杂 


本 文 在 理论 下 


究 时 还 存在 一 些 不 足 之 处 。 首 先 在 静态 分 析 


时 , 只 是 考虑 投资 后 总 成 本 小 于 等 于 投资 前 总 成 本 时 的 AGYV 蔡 


[三 规模 变 大 时 的 AGV 蔡 换 比 。 同 时 ,本文 没 


例 ， 在 以 后 的 研究 中 有 必要 结合 具体 实例 和 
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相关 财务 知识 如 现金 流量 、 行 业 基 准 率 去 考虑 AGV 替换 比 。 
最 后 ， 本 文 只 研究 了 搬运 一 项 操作 的 人 机 蔡 换 ， 在 以 后 的 研究 
中 有 必要 扩展 到 每 一 个 岗位 进行 人 机 蔡 换 研究 。 
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